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ABSTRACT 

SDSS J120602. 09+5 14229.5 is a gravitational lens system formed by a group of galaxies 
at redshift zfg = 0.422 lensing a bright background galaxy at redshift zbg = 2.001. The 
main peculiarity of this system is the presence of a luminous satellite near the Einstein ra- 
dius, that sHghtly deforms the giant arc. This makes SDSS J120602.09+5 14229.5 the ideal 
system to test our grid-based Bayesian lens modelling method, designed to detect galactic 
satellites independently from their mass-to-light ratio, and to measure the mass of this dwarf 
galaxy despite its high redshift. We model the main lensing potential with a composite an- 
alytical density profile consisting of a single power-law for the group dominant galaxy, and 
two singular isothermal spheres for the other two group members. Thanks to the pixelized 



source and potential reconstruction technique of Vegetti & Koopmans ( 2009a ) we are able 
to detect the luminous satellite as a local positive surface density correction to the overall 
smooth potential. Assuming a truncated Pseudo-Jaff'e density profile, the satellite has a mass 
Msub = (2.75 ± 0.04) X lO^^Mo inside its tidal radius of = 0.68'' . This result is robust against 
changes in the lens model, with a fractional change in the substructure mass from one model to 
the other of 0. 1 percent. We determine for the satellite a luminosity of = (1 .6+0.8) x 10^ L©, 
leading to a total mass-to-light ratio within the tidal radius of (M/L)b = (17.2 + ^.5)Mq/Lq. 
The central galaxy has a sub-isothermal density profile as in general is expected for group 
members. From the SDSS spectrum we derive for the central galaxy a velocity dispersion of 
c^kinem = 380 + 60kms"^ within the SDSS aperture of diameter 3''. The logarithmic density 
slope of 7 = 1.7^Q3Q (68% CL), derived from this measurement, is consistent within 1-cr with 
the density slope of the dominant lens galaxy y 1.6 determined from the lens model. This 
paper shows how powerful pixelized lensing techniques are in detecting and constraining the 
properties of dwarf satellites at high redshift. 
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1 INTRODUCTION 

Comparison between numerical CDM simulations and direct obser- 
vations of the Milky Way and Andromeda has shown the existence 
of a strong discrepancy in which the abundance of predicted sub- 
haloes outnumbers that of observed dwarf galaxies (e.g. Kauffmann 
let al.|1993[|Klypin et al.|1999l|Springel et al.|2008[[^avtsov,2010t 
and references therein). Reconciling the luminosity function with 
the mass function is therefore a crucial test for CDM models. With 
the specific aim of addressing this issue, a grid-based Bayesian lens 
modelling code was developed by [Vegetti & Koopmans] ( [2009a| ). 
This technique is able to identify possible mass substructure in the 
lensing potential, by reconstructing the surface brightness distribu- 
tion of lensed arcs and Einstein rings. Several tests on mock data 
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have shown that we can detect mass substructure as massive as 



Msub > 10^ Mo ( Veget ti & Koopmans |2009a|bt 



In a recent application to the lens system SDSSJ0946+1006 from 
the Sloan Lens ACS Survey (SLACS, [Bolton et ([2006}), the 
method has proved to be successful in recovering the smooth 
lensing potential and in identifying a satellite with a high mass- 
to-Ught ratio, (M/L)v > 120 (M/L)v,o, and with mass M^ub ~ 
(3.51 + 0.15) X 10^ Mq, while reconstructing the data to the noise 
level ( [Vegetti et al.||2009| ). However, the complexity of the data 
and systematic eff'ects related for example to sub-pixel structure, 
PSF modelling and spatially varying noise can complicate the 
source and the potential reconstruction and all their eff'ects al- 
ways have to be carefully assessed and quantified. A definitive 
test is therefore required to calibrate the capability of our tech- 
nique. SDSS J120602.09+514229.5 ( [Lin et al.|2 009') with a lumi- 
nous satellite right on the lensed images (see fig. 1) is an ideal sys- 
tem to accomplish this task. In this paper we present a full anal- 
ysis of SDSS J120602.09+5 14229.5, measuring the mass and the 
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mass-to-light-ratio of the dwarf satelUte and showing the strength 
of the method on known satellites. The layout of the paper is as fol- 
lows. In Section 2 we introduce the data. In Section 3 we provide 
a short description of the modelling method and a detailed descrip- 
tion of the main results and in Section 4 we summarise our main 
results. Throughout the paper we assume the following cosmology 
Ho = 73 kms"^ Mpc'^ = 0.25 and Qa = 0.75. 



2 THE DATA 

SDSS J120602.09+5 14229.5 (the Clone) was observed with the 
Hubble Space Telescope (HST) in cycle 16 (P.I.: S. Allam). WFPC2 
images were obtained through three filters, F450W, F606W, and 
F814W. We base our lens model on the F606W data because they 
provide the best combination of depth and resolution. 
Four dithered exposures were obtained, each with an integration 
time of 1 100 s. We retrieved the calibrated exposures from the HST 
archive and used multidrizzle to combine them. The final image has 
a pixel scale of 0.05 arcsec. 

To model the observed structure of the lensed source, we require 
knowledge of the point spread function (PSF) of the drizzled im- 
age. Since there are no suitable stars in the field, we rely on a model 
PSF created with Tiny Tim, v6.3 ( Krist 1993 I^The PSF is gener- 
ated for the position of the lens system on chip 2 of the WFPC2, 
subsampled by a factor 10. Subsampling is necessary because the 
dither pattern involves half-pixel shifts which cannot be taken into 
account by Tiny Tim alone. Instead we rebin the highly subsam- 
pled PSF to the original WFPC2 pixel scale once for each science 
exposure, with the output grid shifted by five subsampled pixels to 
account for half-pixel shifts. The rebinned PSF is inserted at the 
position of the lens galaxy Gl in copies of the four science files 
after setting the image data to zero. These PSF exposures are then 
drizzled with the same parameters as before to create an approxi- 
mation to the PSF of the drizzled science image. 
Light from the outer parts of the main lens galaxies Gl, G2 and G3 
contributes a few percent to the light at the location of the lensed 
arc images. We subtract de Vaucouleur models of these galaxies, 
generated with GALFiij^ version 3.0 (Peng e t al.|2002| ). The eff'ec- 
tive radius of Gl is determined at 3.9+0.1 arcsec. 
Finally, we need to remove the satellite galaxy G4. Since this 
galaxy sits on top of an arc image and it is not a priori clear what 
the background level due to light from the arc is, no unambiguous 
model of G4 can be determined. Due to its compactness we opt 
to subtract a simple Gaussian model with FWHM 0.173 arcsec and 
normalisation determined from visual impression of the image after 
subtraction. Changes in the normalisation within a plausible range 
do not significantly aff'ect the potential reconstruction. 
A spectrum of Gl is available from the Sloan Digital Sky Survey. 
This spectrum has a fairly low signal-to-noise ratio of S IN ^ 8 
(per pixel) and no kinematic measurements are given in the SDSS 
database. With the template fitting method described in Czoske et 
al. (2010, in preparation; see also Czoske et al. 2008) and the Indo- 
US spectrum ( |Valdes et al.|2004| ) of the K2III star HD 195506 as 
template we obtain a good fit (reduced = 1.11; Fig. [2]) with a ve- 
locity dispersion of cTkmem = 380 + 60 km s~^ within the SDSS aper- 
ture of diameter y . This is lower but not inconsistent with the value 
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Figure 1. Overview of the lens system SDSS J120602.09+5 14229.5. This 
false-colour image was created from HST/WFPC2 images in F450W, 
F606W and F814W. 



that'L in et al.| ( |2009"| ) obtained from fitting a singular isothermal el- 
lipsoid model to the lens configuration, criens = 440 + 7kms~^ 
Note, however, that this value applies to the entire mass doing the 
tensing, whereas our value is for the main lens galaxy Gl only. 



3 LENS MODELLING 

The lens modelling is performed using the Bayesian adaptive 
method of Vege tti & Koop mans^ ('2009a') to which we refer for a 
detailed description. Briefly we proceed as follows: 

(i) Initially, we only assume a smooth analytic model for the 
main lens potential (i.e. we only consider Gl, G2 and G3) and 
maximize the relative posterior probability in terms of their lens 
parameters. At this point we ignore the satellite G4. 

(ii) We fix the lens potential at the maximum posterior values 
found in the previous iteration and run a simultaneous pixelized re- 
construction of the source surface brightness distribution s and the 
potential correction 6^. This leads to the detection and localisation 
of possibly present mass substructures in the lens potential. 

(iii) Finally, we build a composite analytic model in which both 
the main lenses and the detected satellite have a power-law (PL) 
density profile and we optimize the relative penalty function for the 
corresponding parameters. The power-law is truncated in the case 
of the satellite. 

In the next sections we describe this procedure in more detail as 
applied to the lens SDSS J120602.09+5 14229.5, and show that it 
is able to detect and quantify the satellite G4. 

3.1 Smooth potential parametric reconstruction 

We initially start with a smooth model that explicitly excludes the 
satellite G4. We model the tensing potential as the combination of 
a single power-law ellipsoid for Gl and two singular isothermal 
spheres (SIS) for G2 and G3 with a surface density in terms of 
critical density ( [Kormann et al.|1994] ) 
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Figure 2. SDSS spectrum of galaxy Gl is shown in the upper panel. The spectrum of K2III star HD 195506 is convolved with a Gaussian line-of-sight velocity 
distribution of dispersion 380kms~^ overlaid in red. The residuals of the fit are plotted in the lower panel, with the expected noise spectrum overlaid in blue. 
The masked regions are for B aimer lines, metal lines that typically show abundance anomalies compared to galactic template spectra, atmospheric absorption 
features and regions that appear contaminated with strong spikes. 



where r = 

Given the relatively high dynamic range of the lensed image surface 
brightness distribution, the source is reconstructed on a Delaunay 
tesselation grid that is built from the image plane by casting every 
second pixel in RA and DEC back to the source plane. The area 
of each triangle (i.e. the grid resolution) depends on the local lens 
magnification. A curvature source regularisation is adopted ( |Vegetti| 
|& Koopmans|2009a| . The free parameters for the posterior proba- 
bility maximization are the lens strength b, the position angle 6, the 
axis ratio q and the density slope y for Gl, the lens strength for G2 
and G3, the strength of the external shear Fsh and its position angle 
^sh. The best PL+2SIS model is reported in Table 1. We find that 
this is an incomplete and simplified description for the true lensing 
potential of SDSS J120602.09+5 14229.5 and do not expect it to 
provide a good description of the data. It does however provide a 
sufficient starting point for the next modelling step. Following |Lin| 
|et al.| ( [2009] ) we also tried a simplified model where the lensing po- 
tential of Gl, G2 and G3 is described by a single global SIE (plus 
external shear), with free parameters being 6, q, the centroid co- 
ordinates Xc and ijc, Fsh and ^sh (see Table 1). For this model our 
results are consistent with [Lin et al.| ( |2009] ) but still provide an ap- 
proximate description of the lens data. Assuming the total mass in- 
side the Einstein radius and the SDSS luminosity-weighted stellar 
velocity dispersion as two independent constraints, we determine 
an effective logarithmic density slope, based on spherical Jeans 



modelling (see Koopmans et al. 2006 for details), of y = l-^-o so 
(68% CL), assuming orbital isotropy (/? = 0), an effective radius of 
3.9+0.1 arcsec for Gl and a seeing of 1.5 arcsec. Conversely, the 
best PL density slope of Gl, i.e. y - 1.58, predicts a dispersion of 
340 km/s. Hence the best PL model is in excellent agreement with 
the measured stellar velocity dispersion, although the SIE model, 
which predicts a large dispersion of 450 km/s is still marginally in 
agreement given the larger error on its measured value. 



3.2 Satellite detection 

The next step is to test whether the pixelized technique is able to 
identify the satellite G4. We fix the lens parameters of Gl, G2 
and G3 to the values found in the previous section and run a lin- 
earized grid-based reconstruction for the source surface brightness 
and lens potential corrections. To ensure the linearity of the solu- 
tion, both the source and the potential corrections are initially over- 
regularized and then the relative regularization constants are slowly 
lowered. Curvature regularization is used for both the source s and 
potential corrections 6^. 

The potential corrections are reconstructed on a regular Cartesian 
grid with 81x81 pixels and a pixel scale of 0.12''. Via the Poisson 
equation 6^ can be translated into convergence (surface density) 
corrections 6k - \{6\l/ii + 6\l/2i)- A strong positive convergence 
correction is found at the exact position of G4 (see Fig.|3]). Smooth 
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Table 1. Best recovered parameters for the mass model distribution for the lens SDSS J120602. 09+5 14229.5. For each of the considered models we report 
the best recovered set of non-linear parameters and the galaxy centroids. 



Model 


Lens 




e 

(deg) 


q 


(arcsec) 


(arcsec) 


r 




6'sh 
(deg) 


logX 


PL+2SIS 


Gl 


1.59 


-50.7 


0.75 


-0.54 


-0.13 


1.48 


-0.06 


-27.8 


48956.90 




G2 


0.54 




= 1.00 


0.38 


-0.98 


= 2.0 










G3 


0.55 




= 1.00 


-2.23 


0.33 


= 2.0 








SIE 


G1+G2+G3 


3.67 


-71.2 


0.76 


-0.45 


0.002 


= 2.0 


0.005 


-63.1 


50267.43 


PL+2SIS+PJ 


Gl 


2.27 


-79.2 


0.80 


-0.54 


-0.13 


1.58 


0.03 


3.53 


102308.08 




G2 


0.17 




= 1.00 


0.38 


-0.98 


= 2.0 










G3 


0.13 




= 1.00 


-2.23 


0.33 


= 2.0 










G4 


2.78 X IQ^^Mq 






3.12 


-2.10 










SIE+PJ 


G1+G2+G3 


3.78 


-74.7 


0.80 


-0.45 


0.002 


= 2.0 


0.02 


-74.1 


114035.73 




G4 


2.75 X \{)^^Mq 






3.12 


-2.10 











non-negligible density corrections are also found on the upper side 
of the arc. These could be related to the fact that the source regu- 
larization is not at its optimal level, but slightly over-regularized, 
or more likely to smooth deviations of the starting model from the 
true mass distribution (see e.g. Barnabe et al. 2009 ). 
Reconstructions with different values of source and potential reg- 
ularization lead to very similar results (see fig. [4|. As expected, 
the potential correction at the position of G4 becomes more ex- 
tended and less concentrated for higher levels of regularizations 
(^ = 3.0 X 10^ and X^^ = 3.0 x 10^ = 3.0 x 10^ and 
X§^]J - 3.0 X 10^) but is otherwise the same for all other combi- 
nations of these parameters. This indicates the robustness of the re- 
sults against changes in the source structure and potential smooth- 
ness. 

Also the single global SIE leads to a similar convergence map, with 
the density correction corresponding to G4 located at the same 
position and having a comparable intensity as for the multiple- 
component model. The satellite detection is therefore robust against 
diff'erent choices for the initial smooth global lens potential. In fact 
a SIE-hPJ model is slightly better than a PL-h2SIS+PJ one. This 
could be interpreted as due to the presence of a common halo for 
this group of galaxies. It is important to note that the convergence 
correction is located exactly at the position of the peak of the sur- 
face brightness distribution of G4 as recovered in Section 2 via a 
Gaussian fit. 

We conclude that the extra freedom allowed to the lens potential via 
the linear potential corrections compensates/corrects for the inade- 
quacies of the global lens potential and both identifies and precisely 
locates possible mass substructure. 



3.3 Satellite mass 

In this section we further quantify the pixelized substructure by an 
analytic model and constrain the relative parameters in the context 
of that model. We assume an analytic mass model consisting of a 
single PL for Gl, two SIS for G2, G3 and a Pseudo-Jaff^e (PJ) for 
G4 as well as a simplified model SIE-hPJ. The Pseudo-Jaff'e profile 
reads as ( |Dalal & Kochanek|2002[|Vegetti et aLpjOQ, ): 

E(r)=^[r-^-(r^ + r?)-^/2], (2) 



where = is the tidal radius for a lens strength ^sub- The 

satellite G4 is centred on the position where the peak of the con- 
vergence correction was found by the pixel-based reconstruction. 
The free parameters for Gl, G2, and G3 are the same as before, 
while the only free parameter for G4 is the mass within the tidal 
radius Msub = ^^^^sub^^c- The recovered best parameters are listed 
in Table 1 for both models. The inferred substructure mass is not 
strongly aff'ected by small changes in the substructure position; a 
systematic change of 1 pixel in the centre coordinates leads, for 
example, to a change in the substructure mass of only 1 percent. 
They respectively lead to a satellite mass and tidal radius 
Msub = (2.78 + 0.04) X lO^^Mo, r, = 0.68" (PL-h2SIS+PJ) and 
Msub = (2.75 + 0.04) X lO^^Mo, r, = 0.81" (SIE-hPJ). 
The reader should not be tempted to compare the diff'erent models 
in terms of the Likelihood reported in Table 1; models can only 
be compared in terms of the Bayesian evidence, which requires 
to integrate over the multidimensional space of the posterior 
probability density distribution over the free-parameters. The 
model comparison is not relevant for our current analysis and this 
step is therefore not carried out. 



3.4 Satellite mass-to-light ratio 

Finally, we estimate the luminosity of G4 by integrating the Gaus- 
sian model to the F606W surface brightness profile obtained in 
Sect. 2. We expect this to lead to a underestimate of the luminos- 
ity, because of the sharply dropping wings of the Gaussian model. 
Fitting more realistic models is, however, difficult due to the com- 
pactness of G4 and the contamination with the arc light. The colour 
of G4 is consistent with that of the main lens galaxies Gl, G2 and 
G3, which indicates an old stellar population. The absolute rest- 
frame B-band magnitude is obtained following the prescription of 
|Treu et aL] { [1999) ) for an elliptical galaxy and is Mb = -17.5, cor- 
responding to a luminosity of = (1.6 + 0.8) x 10^ L©. The large 
error estimate includes the uncertainty due to arc light contaminat- 
ing G4 and the model uncertainty for the light profile. The total 
mass-to-light ratio of G4, inside the tidal radius, is thus (M/L)b = 
(17.2 + S.5)Mq/Lq. As explained, this should be really only con- 
sidered an upper limit to the true mass-to-light ratio. Plausible de 
Vaucouleur profiles are typically 0.8 mag brighter than the Gaus- 
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Figure 3. Results of the pixelized reconstruction of the source and lens potential corrections. The top-left panel shows the original lens data, the middle one 
shows final reconstruction while the top-right one shows the image residuals. On the second row the source reconstruction (left), the potential correction 
(middle) and the potential correction convergence (right) are shown. 



sian model, leading to a total luminosity of Lb = (3.3 + 1 .6) x 10^ L© 
and a mass-to-light ratio of (M/L)b = (8.2 + 2.6)Mq/Lq. This re- 
sult is consistent with little to no dark matter inside the tidal radius 
of this satellite; this is also in agreement with the typical stellar 
mass-to-light ratio at this redshift (M/L)b ~ 5Mq/Lq ( |Treu et aT] 
12005] ). 



4 SUMMARY 

We have applied the grid-based Bayesian lensing code 
by [Vegetti & Koopmansj ( |2009a| ) to the lens system 
SDSS J120602.09+5 14229.5, which has a known luminous 
satellite located on the lensed arc. We have shown that the pertur- 
bation of the lensed arc, created by the satellite, can be used to 
gravitationally identify the satellite itself and determine its lensing 
properties, in particular to get an accurate mass measurement. We 
performed several tests that show that the satellite detection and its 
recovered mass are robust against changes in the source structure, 
level of lens potential smoothness and choice of the smooth global 
lensing model. The main results of this work can be summarised 
as follows: 

• A relatively complex model, containing one single power-law, 
two singular isothermal spheres and a Pseudo-Jaffe satellite, yields 
a satellite mass Msub = (2.75+0.04)xlO^^Mo inside the tidal radius. 
This result is consistent with a simpler SIE-hPJ model. 



• The satellite has a total mass-to-light ratio within the tidal ra- 
dius of (M/L)b ~ S.OMq/Lq, consistent with the presence of little 
to no dark matter inside the tidal radius, assuming a typical stellar 
(M/Lb)^ ^ 5.OM0/L0 

• Gl, the main galaxy in the group, has a density profile which 
is sub-isothermal with slope y = 1.58 + 0.1. This is not unexpected 
for galaxies in groups (e.g jSand et al.|2004] ) 

• We measure for Gl a velocity dispersion of cTkmem = 380 + 
60kms~^ within the SDSS aperture of diameter This is consis- 
tent with the (Tsm value from Li n et al. (2009J obtained by fitting a 
singular isothermal ellipsoid model to the lens configuration. From 
a more proper lensing and dynamics model we predict a stellar ve- 
locity dispersion of 340 km/s for the best PL model of Gl that as 
a logarithmic density slope of 7 = 1.58. Conversely, we predict 
a density slope of 7 = ^-^-ofo (68% CL) from the observed stel- 
lar velocity dispersion. This agrees very well with that determined 
from the PL model of Gl, but is also still marginally in agreement 
with the SIE model. 



This paper demonstrates the great potential of pixelized lens- 
ing techniques in robustly identifying and measuring the key prop- 
erties of small mass structure/dwarf satellites in distant galaxies. 
The application of this method to a large uniform set of lens galax- 
ies will allow in the near future to constrain the general properties 
of mass substructure in galaxies and to test the CDM paradigm on 
these small scales. 



© 2002 RAS, MNRAS 000,[l]-?? 



6 S. Vegetti, O. Czoske & L. VE. Koopmans 





\ 






-4 -2 

arcsec 


2 4 




^ 






-A -2 

arcsec 


2 4 


\ 









I 











-4 -2 


2 4 
arcsec 








f 






-4 -2 


2 4 
arcsec 




j 







I 











-4 -2 


2 4 
arcsec 








V 




T 

Y 






-4 -2 


2 4 
arcsec 



I 



I 




Figure 4. Results of the pixelized reconstruction of the convergence corrections for different values of the source and potential regularization Xs - 3x10^ (top 
row), = 3 X 10"^ (middle row), = 3 x 10^ (bottom row) = 3 x lO'^ (left column) and /l^^ = 3 x 10^ (middle column) and = 3 x 10^ (right column). 
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